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We extend Gor'kov theory to address superconducting pairing at high magnetic fields and gen- 
eral temperatures with arbitrary attractive interaction strength. This analysis begins with a new 
interpretation of the high-field Gor'kov gap equation which we associate with an instability in a 
generalized particle-particle ladder series. Importantly, this interpretation of the non-linear gap 
equation enables a treatment of pairing which is distinct from condensation. We also show how to 
consolidate two distinct fermionic pairing schemes in real and momentum space, both corresponding 
to an Abrikosov lattice. Numerical results for the fermionic local density of states demonstrate that 
gapless structure in a field is robust and presumably relevant to quantum oscillation experiments. 
We find that despite their differences, both pairing schemes contain very similar physics. Our for- 
malism is designed to explore a variety of magnetic field effects in the so-called pseudogap phase 
and throughout the BCS-BEC crossover. 



I. INTRODUCTION 



High magnetic field superconductivity is a difficult 
theoretical problem with important implications for ex- 
periment. Observations of quantum oscillations in the 
superconducting phase of conventional superconductors 
were initially unexpected^ and further unusual oscilla- 
tions have been observed in undcrdoped cupratesP At 
the same time, theoretical treatments reached surprising 
conclusions, with some investigations finding re-entrant 
superconductivity at extremely high magnetic fields. 3 
Adding to the complexity is the observation that the in- 
troduction of Landau levels into a pairing scheme appears 
to cause a three-dimensional superconductor to behave 
like a one-dimensional system. This greatly enhances 
fluctuation effects^ and may even de-stabilize superfiuid 
condensation. 

The goal of the present paper is to set up a founda- 
tion for addressing these issues by extending the standard 
Gor'kov (BCS-based) approach to apply to general tem- 
peratures T, away from the instability regime, and to 
stronger-than-BCS attractive interactions. Although the 
immediate focus of this paper is on the ordered phase, we 
use the Gor'kov theory extension to arrive at a compati- 
ble description of the normal phase as well. Importantly, 
this normal phase may possess a rich structure associated 
with precursor pairing in the presence of magnetic fields. 

In a related pape^ we focus on this disordered phase 
and explore the notion that some degree of pairing at the 
onset of condensation may be necessary to avoid a strict 
dimensional reduction that prohibits condensation into 
a superconducting state altogether. Furthermore, these 
excited pair states may be present in systems such as high 
temperature superconductors or in fermionic gases in the 
BCS-BEC crossover regime, where the non-condensed 
pairs are associated with a pseudogap state. It should 



be stressed that the non-condensed pairs we consider are 
distinct from conventional fluctuations. These pairs arise 
from strong attractive interactions, not from low dimen- 
sionality (and/or disorder effects) which give rise to con- 
ventional fluctuations. Our BCS-BEC based approach is 
similarly distinct from the so-called "phase fluctuation" 
scenario which is based on soft phase fluctuations pre- 
sumably arising from low plasma frequencies. Indeed, 
since we are contemplating both charged and uncharged 
superfluids, the issue of soft plasma frequencies is not 
particularly generic. 

Gor'kov theory addresses the fermionic degrees of free- 
dom. Two proposals^ have been put forth to describe 
the nature of those fermionic pairs which form the con- 
densate in the presence of high magnetic fields. These are 
associated with orbit center-based and magnetic lattice- 
based pairing schemes. Here we show how the physical 
implications of each are similar and that both lead to gap- 
less fermionic states which are thought to be the basis for 
observed quantum oscillations!^ A central contribution of 
the present paper is to demonstrate that the (analytically 
tractable) theoretical structure of these different pairing 
approaches can be consolidated into a more general for- 
mulation which addresses the non-linear structure of the 
Gor'kov theory. This is in contrast to a substantial frac- 
tion of the literature on high magnetic field supercon- 
ductivity which deals with the linear regime where the 
gap is small. We thus arrive at an interpretation of the 
Gor'kov gap equation which allows us to extract a set of 
particle-particle ladder diagrams which properly charac- 
terize the pairing fluctuations or non-condensed pairs in 
the presence of high magnetic fields. 

We begin in Section [TT] by deriving the Landau level 
representation of the Gor'kov equations. In Section III 
we discuss a "diagonal" approximation made to these 
equations which is thoughtP to be suitable for high mag- 
netic fields. Following this, in section [TV] we focus on the 
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gap equation and discuss its relationship to a divergent 
particle-particle ladder series. In section [V] we show that 
this divergent series captures the contribution of non- 
condensed pairs to the theory. 

The second half of the paper is less general and more 
concrete. Here we focus on implementing and compar- 
ing different state-space representations of this pairing 
theory. In Section |VI| we discuss the two different ex- 
isting implementations. In Section |VII| we show how to 
use a "tight-binding" approximation to orbit-center pair- 
ing to make the theory analytically tractable, facilitating 
comparison of the two implementations. Finally, in Sec- 
tion [Vlll] we show how gaplessness is robust among these 
pairing theories by comparing local density of states cal- 
culations. Our conclusions are presented in Section |IX| 



II. DERIVATION OF GOR'KOV'S EQUATIONS 
IN A LANDAU LEVEL BASIS 



We begin with the Gor'kov equations formulated in 
real space in terms of the gap A(r) and Green's functions 
in integral form: 



G(r, r'; iu) = G (r, r'; iu) - J dr"dr'"G Q (r, r";iu) 
x A(r")Go(r"',r";-iu)A^(r"')G(r"',r';iu) (1) 

At(r) = -|^ / dr'G(r',r;iu)G (r',r;-iw)A^(r') 

iuj 

where g is the interaction strength, and iu (iQ) denote 
fermionic (bosonic) Matsubara frequencies with the tra- 
ditional subscripts omitted for clarity. These equations 
allow the identification of the self-energy E, 

Y,{r,r'-iu) = -A(r)A^r')G {r',r:-iu), (3) 



and we now rewrite Eqs. (|l|3j) in terms of the Landau 
level representation that diagonalizes the non-interacting 
Hamiltonian "Ho^The bare Green's function Go is then 
given by 



n 



iu - £„ 



(4) 



where n = (N,p,k z ) is the Landau level state, with TV 
the Landau level, p the degenerate Landau level index, k z 
the momentum in the z-direction (parallel to a magnetic 
field B), and the energy of a particle in state n. G, 
however, is not in general diagonal in the Landau level 
representation, and is given by 

G(r,r'-iu) = ]T G mm ,{iu)i> m {r)^ m ,{r'). (5) 



In this representation, Eq. ([I I is multiplied by 
i } m{ r )'4 } rn'{i'') and integrated over both r and r' to give 

G mm ,{iu) = G° m (iu)S mm , -J2 J dr"dr"'G Q m {iu) 

In 

x^(r")A(r")GO(-z W )^(r-"')^ t (r-")A t (r"') 
x G 

nrn (iu)ip n (r"'). 
We then define a "state-space gap" 

A m „ = J drA{r)if;i(r)iPi(r) 

and obtain 

G mm >(iu) = G° m (iu)8 mm , - y^G^(zo;)An 



(6) 



x G¥(-iu)A\ n G nm ,(iu). 



(7) 



We multiply the gap equation, Eq. (J2j|, by A(r), and 
express the right hand side in the Landau level represen- 
tation to find 

Jdr\A(r)\ 2 = -| E E/ dr'drG mm ,(iu) 

mm'n iuj 

xG° R H u )At(r')A(r)fc(r')iWi(r>lW. 

Using the above identity for A mn , the gap equation be- 
comes 

i =- g - a y y Am ' nA ™ n 2 G mm ,(iu)G° n (-iu). (8) 

P ^ ^ dr\A(r)\ 2 V ; nV ' W 

mm'n iu) J I V /I 

The self-energy is similarly expressed in the Landau 
level representation, 

Z(r,r';iu) = ^ ^mm'{i^m{r)^ m ,{r'). (9) 

mm' 

Using Eq. ^ for S, we find 

E w (tw) =-Y J G° n {-uu)A mn Al ml . (10) 

n 

In summary, the Gor'kov equations in the Landau level 
representation for a constant magnetic field are 

G mm >{iu) = G a m {iu)5 mm , - y^G"„(z^)A mi 



In 



x Gf(-iu)Aj n G nmt (iu) 



(7) 



9 \ ^ \ " ^■m'nA\ nn 



(3 ^! 4^ fdr\A(i 

mm n iuj J 1 v 



G mm >(iu)G° n {-iu) 



V mm ,{iu) = -Y,G a n (-iu)A mn Al m „ 4T0) 
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III. INTRA-EIGENSTATE PAIRING 

To make f urther analytical progress, we mu st simplify 
Eqs. ( |7|8|10| ). Thus, as in the literature^ 121121 we assume 
that the pairing involves degenerate eigenstates. This as- 
sumption is justified if we are in a regime where |A| is 
much less than the splitting between the Landau levels 
hw c , where u c = eH/(Kcm) is the cyclotron frequency. 
This will be a good approximation at very high fields 
where the splitting between Landau levels is large, the 
regime we explore in this paper. This approach has been 
analyzed carefully in Ref. [5] which showed that in this 
regime inter-Landau level effects that we neglect are per- 
turbative. 

We also simplify by defining 



A^(C) = 



A ron (C) _ A m „(C) 



Jdr\A(r,Of 



(11) 



Note that we have introduced an important parame- 
ter £ which labels different functional forms of A(r). 
Throughout this paper we will associate the particular 
value £ with the condensate configuration A sc (r,£ )> 
which may be distinct from configurations occupied by 
non-condensed pairs, discussed below. It is convenient to 
define 



4>mm'(0 



VA° (C)A ot ,(C). 



(12) 



In order that all the potential pairing partners n of a 
state 77i are energy degenerate, the index n is in the same 
Landau level as m, and pairing occurs between states 
with z- momenta fc z and —k z . This allows us to write the 



fermionic self-energy Eq. ( 10 ) as 



S w (iw) = ~G° N (k z ;-iu) |A| 2 0w(Co), (13) 

where we write G° in terms of only the Landau level 
N and the z-momentum k z of m, with N m — N m i and 
k Zm = k z , . This simplification of the gap equation, 
Eq. leads to 

1 = "I E J2^ m '(Co)G mm ^)G° N (k z ;-iu>). (14) 



mm' iu) 



IV. THE GAP EQUATION AS DIVERGENT 
PARTICLE-PARTICLE LADDER SERIES 



The above formulation enables a consol idation of stan- 
dard approaches in the literaturePZEEfi] j n this context 



we re-interpret the gap equation of Eq. ( 14 ) as a gener- 



alized Thouless condition, which applies to all tempera- 
tures below the transition. In the process we show that 
this gap equation serves to identify a particular particle- 
particle ladder series. The divergence of this series is 
a prerequisite for, and an indicator of, the supercon- 
ducting state. Importantly this ladder series also leads 



us to a characterization of the associated non-condensed 
pairs, which may form above the critical temperature for 
stronger-than-BCS interactions. 

It is essential first to characterize the degrees of free- 
dom available to these non-condensed pairs. In the z- 
direction parallel to the magnetic field, the system be- 
haves as in zero-field, where condensed pairs are con- 
structed of fermions with momenta k z and — k z . Thus, 
excited pairs must have nonzero total momentum, and we 
can describe the general pairing of non-condensed pairs 



as being between momenta k z and — k z 



rm 



In the plane perpendicular to the magnetic field, the 
condensed electrons are those which pair to form the real- 
space superconducting gap. Here and throughout we dis- 
tinguish the order parameter A sc (r,£o) from the excita- 
tion gap. In a mean- field scheme, such as ours, where 
individual vortex fluctuations are not included, we as- 
sume that (o corresponds to an Abrikosov lattice with 
functional form A°(r,Co) = A sc (r, Co)/A. We take the 
non-condensed pairs to be those that form other real- 
space gaps A°(r, C) for C 7^ Co- Finally, we also allow the 
Matsubara frequencies to appear with total frequency iQ. 

We next introduce the pair susceptibility, 



a y y y y ^mm 



( l > li m '(C)G m m'(iuj)G N (k z - q z ;iQ,-iu)). (15) 



mm' iuj 



With this important definition, and the introduction of 
non-condensed pairs, we return to the gap equation, 
Eq. ( 14 ) , which can be rewritten as 



1 + 9X 



(o,qz=0;in = 



0. 



(16) 



One can interpret this equation as reflecting a divergence 
of a particle-particle ladder summation, shown in Fig. [T] 
We argue below that the gap equation is to be associ- 
ated with a i-matrix formed from the ladder diagrams in 
Fig. [lj given here in an abbreviated form by 



t^(c,q z;i n) = 



9 

9X{C, a.* 



in) 



(17) 



(see Section [V] for details) . For the condensate config- 
uration, £ = Co, q z = Q = 0, the yj-matrix thus di- 
verges for all temperatures below the instability, as in 
a Bose Einstein condensation condition, where the pairs 
have vanishing chemical potential. This ladder diagram 
set is to be distinguished from a series which was previ- 
ously identified to correspond to the specific instability 
point Ref. [T^l Here, the condition is stronger, as the in- 
corporation of one dressed G and one bare G° makes 
this summation valid throughout the superconducting 
regime. Furthermore, as we explain in Section [Vj it can 
be extended to the entire strong attraction (BCS-BEC) 
regime where there is a finite excitation gap at the in- 
stability, and in w hich pairing and condensation must be 
distinguished! 2 ^ 
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FIG. 1. The particle-particle ladder summation corresponding to a f-matrix i pg , the divergence of which gives the general gap 
equation for the superconducting state. 



V. CHARACTERIZING THE 
NON-CONDENSED PAIRS IN THE GOR'KOV 
EQUATIONS 

BCS theory represents a very special case of super- 
fluidity in which pairing and condensation take place at 
the same temperature. We have just argued that the 
Gor'kov gap equation of Eq. ( [l6| and the closely related 
i-marix of Eq. (17) effectively constrain the nature of 



non-condensed pairs which below the transition condense 
into a state with gap structure A°(r, (o). 

For the moment, we consider Eq. (17 1 as an appropri- 



ate characterization of the i-matrix (or effective propaga- 
tor) associated with the non-condensed pairs. We next 
characterize their feedback into the Gor'kov equations, 
by including t in the self energy. In a strict Gor'kov the- 
ory £w(iw) = 

E <f^ m '(0^ m >(<,9z]i^)(^{gx-k,iin-iu>), (18) 



where 



i sc (C, g z ; iQ) = -S(C - Co)%,)*(ifi) |A s f , (19) 



with A sc corresponding to the gap A in Eq. (11 1 



We now include in the self-energy the non-condensed 



pair propagator t pa given by Eq. (17), and in this way 
go beyond strict Gor'kov theory. This contribution is 



<f> 2 mm ,(()t Pg (<;,<lz;M)G N (q z - k z ;in - iw). (20) 



Furthermore, since 1 + <7x(Co, 0; 0) diverges below the 
critical temperature, t ps will be strongly peaked for pa- 
rameters ((,q z ;ifl) « (Co,0;0). Since iQ and q z will be 
small for the primary contributions to i ps , we approxi- 
mate the right hand side using G° N {q z — k z ;iQ — ioj) m 
G%{-k z ;-iw). 

Then the total self-energy is now 

= G° N (-k z --iw) £ 4> 2 mm >{Q 
x (i sc (C,<7,;^)+^(C,g z ;^)), (21) 



where S^ OT , is the original self-energy from Gor'kov the- 
ory defined in Eq. ([To]) , and E£f m , is defined in Eq. ((20 1 
We can further simplify by defining a total |A| 2 , 



lAr = ia s 



l A pg 



with 



i A Pg 



y 9 



(22) 



(23) 



and <#, m ,(C)) c = 

E c , 9iji o (* SC (C, *«) + * pg (C, g z ; *fi)) 



(24) 



This leads to an expression for the self-energy which can 
be written compactly as 

V mm ,{iu>) = -G%(-k z ; -iu) (<fLn'(0) c l A | 2 ■ (25) 

Together with the numb er equation N 



2 J2 m ,iu G mm {iu) and Eqs. |7j8| this forms a sys- 



tern of equations which can be solved self-consistentlyP»Q 
In this way we have modified the Gor'kov theory to 
extend BCS theory into the more general regime where 
pairing and condensation are treated differently. 

Finally, we turn to a more detailed interpretation of 
the ladder diagrams and the related pair susceptibil- 
ity. In Fig. [T] particles forming a pair may interact to 
create a new non-condensed pair in the same excited 
pair state. That is, we consider a pairing propagator 
i™" ' pg ((,q z ;in) between states m,n and m',n', where 
m = (N,p, k z , iuj) and n — (N' ,p' , q z — k z ,ifl — iui) are 
pairing partners for a given real-space configuration £, 
total ^-momentum q Zl and total Matsubara frequency 
iVt, and m! = (N",p" , k' z ,iw) and n' = (N"',p"',q z - 
k' z ,in — iuj') are pairing partners sharing the same £, 
q z and in. Note that m now includes frequency, and 
intra-eigenstate pairing is not enforced here. The bare 
interaction between the two pairs, as an extension of 
the BCS mean-field Hamiltonian, is in a separable form 
gV^ n ' = 5 A£t n (0A^„,(C). Then the infinite particle- 
particle ladder summation is 



m n 
mn 



C„ B,pg (C,g.;in) = sK 

9V™y'G n „G m „ m ^'>™(t, q z ;in). (26) 



5 



As V is separable between m, n and m' 
summation has the solution 



a.vci n ,pg 
mn 



(C,q z ;ifl) = 

gv r 



m n 
mn 



9 Em" 



m" \m'" ,n" m/ n n 



\rm"n" r^O f* 



, this infinite 



(27) 



Using the expression for V, and applying the intra- 
eigenstate pairing approximation, this further simplifies 
to 



j.m n ,pg 



(C,«.;iO) = 



gA% n (QA° m , n ,(Q 

l " m "'{OGm"m'"G° N (k z 



q z ;iCl — iui) 



This in turn gives the i-matrix appearing in the self- 
energy as 



fe; iQ) = <^ ro ,(C)i pg (C,<7*;^) (28) 



with i pg (C, g z ;i^) defined in Eq. (17). Note this also 



justifies our definition of \ m Eq. (15 



In summary beyond the weak attraction limit, pairs 
can form above the superconducting transition temper- 
ature. These pairs of electrons form in excited states 
rather than the quantum ground state. Interpreting the 
Gor'kov gap equation as a Bose Einstein condensation 
condition allows us to specify a precise ladder series struc- 
ture for these non-condensed pairs. Because it will take 
us too far afield, elsewhere, we discuss a precise treatment 
of the parameter £, in more microsopic detailP 



VI. SPECIFYING THE PAIRING BASIS 

We now turn to the specific pairing choices that can 
be made to solve the system, and explore their physical 
effects. We now consider only a single (, and for nota- 
tional simplicity we now omit this parameter. Ensembles 
of £ are discussed elsewhereP In order to make further 
progress, we must specify A(r) as well as the Landau 
level state basis we use. This requires that one deter- 
mine the pairing partners of e ach fe rmionic Landau level 
state. Based on previous worlPSllS here we analyze and 
compare two complementary choices for the case of an 
Abrikosov lattice. Throughout we use the Landau gauge 
A = (0, Bx, 0) in which the energy gap of an Abrikosov 
lattice, A(r) is given by 

A(r) = C^exp (^mA ^., (r) 
m * ' 

V2b x 



L x L y L zy /TTl H 



J^exp 



y 2 i 2imb x y (x - mb x Y 



VK—m + > 
a l 2 H 



I 2 



(29) 



where V'jvx k ^ s the orbit-center Landau level state for a 
charge-2e particle with Landau level N, orbit center X, 
and z-momentum k z . The normalization of A(r) is cho- 
sen such that J dr\A(r)\ 2 — |A| 2 , with L x , L y , and L z 
the sample lengths and In = y/hc/eH the Hall length. 
Finally, the Abrikosov lattice is characterized by unit vec- 
tors a = (0, a, 0) and b — (b x ,b y , 0). Ijj is related to the 
unit vectors by ab x — For a square lattice b = 

(a, 0, 0) while for a triangular lattice b = (yia/2, a/2, 0), 
and in general we capture all Abrikosov lattices using 
the method in Ref. Q21 C = b y /a + ib x /a. Throughout 
the rest of this paper, we implicitly take b x , b y , and a 
to be functions of C through the above formula and the 
restriction that ab x — t^I\- 

The two natural choices are to use either orbit-center 
wavefunctions, in which case orbit centers positioned 
symmetrically about a lattice site pair together, or to 
use magnetic translation group wavefunctions, in which 
case wavefunctions with opposite reciprocal lattice vec- 
tors are paired. We describe those pairs in terms of the 
notation of A Q mn (£) and </> mm /(C) defined above. 



A. Orbit-center Pairing 



One choice of pairing, originally presented by Ryan 
and Rajagopal is to have fermions pa ir ab out lattice 
site positions in the orbit center basis! 6 ! 14 " 16 ! Because 
the bosonic wave functions forming the Abrikosov lat- 
tice are positioned at orbit centers X = mb x , fermions 
which are equally spaced apart from these positions can 
pair, so that the pair wave function can be represented as 
(r) = 



T pair 

* Ni,N 2 ,mb, 



,Y,k, 



tfr iV™ ft" +Y,k z ,-\{ r ) V' JV™™ - Y, - k z , 1 ( r ) • 



Here 



/ fcrmion 



(r) 



L y L z 2 N N\ 



1/4 



x exp ik z z + 



iXy {x-Xy 



I 2 



2/ 2 



H 



N 



X 



(30) 



where Hn is the Nth-order Hermite polynomial. The 
A mn (C) an d associated factors are calculated in Ap- 
pendix [A] Quite generally (presuming inter-eigenstate 
pairs), the result is 
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(a) 
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FIG. 2. (Color online) Schematic version of t he two pairing models, (a) the restricted "tight-binding" pairing for case (i), 
in the case of a square lattice (see Section VII I. The fermionic states (small open red and filled green circles) can pair with 
each other only via the solid blue (closest pairing partner) and dashed purple (next closest pairing partner) lines. These pairs 
in turn connect "nearest neighbor" states via an off-diagonal Green's function matrix element (dotted green lines; equivalent 
lines for the red states are omitted). Furthermore, the pairs can be interpreted as forming bosons at lattice sites (large blue 
circles), (b) The case (ii) version of pairing, in which states (filled green and open red circles) are in a reciprocal lattice of the 
magnetic translation group, and pair with states of opposite momentum index (blue lines). 



(,f j *; a , 2 oJV +JV 1 ^WWe-y/a iiX = mb x , m G Z, 

lX m=(N u X+Y,k z ).n=(N2,X-Y-k z )\^! ~ \ V v , \ SL ) 

otherwise. 



B. Magnetic Translation Group Pairing 



The other choice of pairing, or iginally presented by 
Dukan, Andreev, and Tesanovic ! 7 ! 10 ! and in parallel work 
by Akera, MacDonald, Norman, and Girvinp^EH anc i 
Nicopoulos and KumarJ^Sl is to use the magnetic trans- 
lation group (MTG) for the fermions, with an index 
fcP Following Ref. HQ] we choose a unit lattice site for 
the magnetic translation group, which must be twice 
the size of the Abrikosov lattice site, with unit vec- 
tors 2a and b. Then the reciprocal unit vectors are 



a* = {~b y /l 2 H ,b x /l 2 H ) and b* = (2a/l%,0) so that 
G^a* = 2nSij where are the unit vectors and a* are 
the reciprocal unit vectors. Restricting k — (k x ,k y ) to 
be within the limits of the cell (a*, 6*) gives a complete 
set of functions!^ We also note that ip is now dependent 
on (, the configuration of the Abrikosov lattice. 

This pairing occurs between opposite k. To conform 
with our assumption of diagonal pairing, we need to spec- 
ify that N\ = N 2 , and then 

= V^Et(r)V#™J^(r) 

with 



formion 
N,k,k z 



2 N N\ V nl% 



i/4 



L X LyL; 



exp (ik z z) 



E 



exp 



2a 



v m 2 



x exp 




a / 



2/ 2 









\ Hn { 




- -) H] 




a I J 



/In) 



(32) 



The A^ n (£) and associated factors are calculated in 
Refs.lSlandlPul with the result that for the lowest Landau 



level 



^m=(0,fe,fc r ),n=(0,-fe,-fc»)(0 



L x L v L z \/2t:Ih 



x exp 



(- (kylH) 2 ^j 3 ( \-k x + iky] b x 



—by iitl 2 H 



(33) 
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where 9s(u\t) = Y^n=-oo exp (2inu + mr?T) is the third 
elliptic theta or Jacobi theta function. Further A^ n (£) 
for higher Landau levels can be iteratively calculated 
from thisP^l 



VII. NEAREST NEIGHBOR PAIRING 
APPROXIMATION FOR REAL SPACE PAIRING 

While the MTG method for pairing results in each 
fermion pairing with exactly one other degenerate eigen- 
state, that is not the case for orbit center pairing. If 
X is a lattice site, then a fermion at X + Y can pair 
not only with X — Y, but also with X — Y + 2nb x with 
n G Z. A very complicated matrix thus results for the 
Green's function, but here we demonstrate that this ma- 
trix can be substantially simplified such that it is analyt- 
ically tractable yet incorporates the important physics. 

We begin by noting that A^ n is proportional to 
exp(— Y 2 /l 2 H ) where Y is half of the distance between 
the two states that pair. Thus, pairing between states 
that are well separated will be exponentially suppressed. 
Using only a single pairing partner is inadequate, how- 
ever, because it fails to capture interference effects be- 
tween two pairing partners, as will be demonstrated be- 
low. Using two pairing partners does capture these ef- 
fects and reproduces the correct local density of states for 



the lattice. This is reminiscent of a hierarchical "tight- 
binding" scheme, in that the nearest pairing partner can 
be considered an "on-site" interaction, and the second- 
nearest pairing partner allows for interactions between 
lattice sites. 

To implement this approach, we begin with Eq. [7J 

G mm ' (ibJ) = 

G° m {iu:)5 mml - G° m (icj) A mJ G?(-iw)A| n G nin /(«j). 

In 



This can be further simplified by separately including the 
Green's functions of pairing partners, resulting in a new 
equation for the Green's functions G mm i (iuj) — 



G° m {iuj)8 mn + G° m {iu)) ^ A m „[-Gw(i 



(34) 



Here by G we mean (a) flipping the sign of ioj and (b) 
conjugating the A appearing within G. By this defini- 
tion, we also have that G mm ' (iw) = 

G° m (-iu)5 mm/ + G° m {-iw) ALGw M- (35) 



This correctly reproduces the original Green's function 
above, and by inverting Eq. (34), we have that 



G _1 (iw) = 



A X-Y,X+Y 








A f 
A f 

L ^X+2b x -Y.X+Y 








A X+2b x -Y.X+Y 

-IGx^-yHu)}- 1 

Ax+2fc x -y : x+2b x +y 







A f 

L ^X+2b x -Y.X+2b x +Y 

A f 

^X+4b al -Y,X+2b x +Y 



A 







X+4b x -Y,X+2b x +Y 



'■J 

(36) 



Here we use only the index corresponding to the orbit 
center, which is the only one which varies throughout 
this matrix. We also take Y > 0; if Y < positions in 
the matrix will flip but it will be otherwise unchanged. 
Also note that the values of Gq 1 are actually independent 
of orbit center; the indices remain for clarity. 

This Green's function matrix is computationally sim- 
plified as compared to the full orbit-center pairing 
scheme. Furthermore, in the case of a square or trian- 
gular lattice, its eigenvalues can be found analytically. 
For both a square lattice (£ = i) and a triangular lat- 
tice (C = 1/2 + \/3i/2), we have the property that 
A x-Y+2b^x+Y+2b :c = A x-Y,x+Y> implying that the en- 
tries in G _1 repeat with period two along the diagonal. 
We can then posit a plane-wave solution for the eigen- 



vectors, with a and b constants: 
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e 1 ^ {a\X rn + 2qb x ) + b\X rn + 2qb x )) , (37) 



where for compactness we set fh = X — Y, m = X + Y, 
and n = X — Y + 2b x , using the symmetry properties of 
A° „ . We obtain 



G~ x v = 



\ 



At m a + (w + £)e ifc *6 + A^e^a 



V 



(38) 



/ 



The energy eigenvalues of this system will be 
zero eigenvalues of G -1 . For a nontrivial solution, 




FIG. 3. (Color online) Plots of the local density of states 
N(r;uj) vs. a; on a square lattice, for A sc = 2 and S = 0.01, 
and normalized to a maximum LDOS of 1.0. (Main figure) 
Case (i), using the nearest neighbor approach for calculations 
with L x — 4b x . (Inset) Case (ii). 



we obtain -A' nm (A„ m e 
A*f ( A - _i_ A- p ik i 



£ 2 + IA- I 2 



+A fim ) + - e) - 

lk,\ _ Q_ naV g _ 

2 



A- l z + IA- \*-\-e~ ik iA- A^ + p ,i )At A- 



(39) 

Here A^ m = C, fim e i ^/< 1 )( x / fc *) 2 , where C„ m is a real 
number, and A fim = C fim e i7r ^/ a W x+1 )/ b -] 2 . We further 
simplify 



p-ikjA /vt 

c "mm'-'nm 



o*jAt A- 



2Rp(p _i ^A- At 

1 ° L -*mm L - 1 nm 



2G mm CV> m R,e ( 6 



where we have used the b x and &j, specific to a square or 
triangular lattice, and the + is for X/b x an odd integer, 
— for X/b x an even integer. We then have — uj 2 



VIII. RESULTS AND DISCUSSION 

This paper has been rather extensively devoted to the- 
oretical formalism which characterizes the fermionic de- 
grees of freedom. Among the most direct experimentally 



(a) 



0.5 
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BB 

-1 -0.5 0.5 



0.95 



■ 0.9 

U0.85 
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0.75 



-1 -0.5 0.5 



FIG. 4. (Color online) Plots of the local density of states 
N(r;uj) for (a) case (i) and (b) case (ii), normalized to the 
local density of states at the vortex cores (located at (x, y) = 
(±0.5, ±0.5)) and as a function of position on a square lattice, 
for uj = 0.4, A sc = 2, and 8 = 0.05 (slightly higher than Fig. [3] 
to permit faster computation without a substantial change in 
accuracy) . 



relevant consequences is the behavior of the local den- 
sity of states (LDOS) N(r;co) vs. uj. Here we address 
this density of states in the very low temperature regime. 
This is experimentally accessible using scanning tunnel- 
ing microscopy.^ The local density of states N(r;u>) is 
calculated via N(r;uj) = 2 Im G ret (r, r; uj). We deter- 
mine G^(r, r'; uj) = £„„„, G mjn .{u ± i8)i> m {r)^ m ,{r') 
in the limit 5 — > + , and for simplicity we address the 
square lattice at the lowest temperature where there is 
only a condensate. Also for simplicity our illustrative cal- 
culations are for the limiting case of the lowest Landau 
level, an s-wave gap, N x = 2 for orbit center pairing, and 
normalized to set the mass m = 1. 

It has been argued quite generally that N(r;uj) ex- 
hibits a gapless behavior, in contrast to superconduc- 
tivity without a magnetic field. This gaplessness is 
dueP to the fact that all fermions are delocalized, un- 
like in the vortex cores of the low field limit. This 
observation has direct application to magnetic oscilla- 
tion measurements, as a gapped state would dampen 
these oscillations significantly.^ More specifically in the 
MTG pairing scheme, systematic studies in Ref . ITTil show 



? + \A mm \ 2 + \A nm \ 2 + 2Cf hm C nm Re ( e -**, e ±**<6«/«)) . | A 

(40) 

We also note that for a system of L x = 2N x b x , we must 
have N x kj = 2nj to enforce periodic boundary condi- 
tions, giving kj — 2irj/N x with j ranging from to 
N x — 1. With this approximation in place, we may turn 
to the results of these pairing theories. 



that the energy eigenvalues satisfy Ej^ t k,k z = k + 

Because A mn always fea- 
tures a zero for the Abrikosov lattice, E will as well. 

The situation for orbit center pairing is more complex, 
but with the nearest neighbor pairing scheme above, we 
are able to analytically demonstrate gaplessness in this 
case as well as for the square and triangular Abrikosov 
lattices. Such gaplessness will occur when Y = 0.5b x , 



which means the magnitudes of |A J; 



A* 



and 



A nm A mm 



vious section), 
also have that 



are equal (using the notation from the pre- 
Thus, for gaplessness to occur, we must 



arg 



-ikj ±z7r(& y /a) 



(41) 
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or jgapicss = ^fli^.J For a square lattice, any 

N x — 2m, m £ N will thus feature a gapless state, 
while for a triangular lattice, any N x = 4m will fea- 
ture a gapless state. This also shows the importance of 
including "nearest-neighbor" pairs; neglecting them can- 
not produce gapless states. 

Beyond this analytical assessment, we observe addi- 
tional similarities and differences between the pairing 
schemes, and gain intuition about experiment. Figure 
[3] presents a plot of the local density of states as a func- 
tion of energy for the two different cases and for two 
different positions of the probe: one at the vortex cen- 
ter (dashed) and the other at the point of maximum | A| 
(solid curve). Again prominent in the features of these 
LDOS plots is that both cases show gapless behavior. In 
both cases, there is only one "nodal" state m per lattice 
site at which its total excitation E = 0. This produces 
the distinctive parabola-like shape in Fig. [3J which still 
touches N(r; 0) = but does not exhibit a complete gap 
away from uj = 0. The variation between the solid and 
dashed curves is rather small, also reflecting this point. 

Figure [4] presents a contour plot of N(r;uj) as a func- 
tion of r. This figure nicely illustrates that the real space 
and reciprocal lattice space pairing schemes are rather 
similar here. Both reflect the symmetry of the Abrikosov 
lattice. However the former shows slightly more contrast 
than the latterP^l 



IX. CONCLUSIONS 

It is anticipated that the formalism in this paper will 
be relevant to both ultracold gases in the BCS-BEC 
crossover regime, and possibly to the pseudogap phase 
in high temperature superconductors. Magnetic field ef- 
fects in the latter have revealed a number of mysteries, 
which appear to be associated with a normal state pairing 
gap. For the cold Fermi gases there is considerable inter- 
est in effects arising from artificial gauge fields or rapid 
rotation. Previous work on these Fermi gaseJ^ has pre- 
sumed, rather improbably, that even in the BEC regime, 
pairing and condensation appear at the same tempera- 
ture. Also important is a better understanding of the 
normal state and of how condensation can take place if 
the superconductor or superfluid is effectively one dimen- 
sional. 

This paper has focused on the nonlinear gap region 
in the presence of a magnetic field. A next step is to 
address calculations of the onset of superfluid coher- 



ence at temperature T C (H), which is taken to be less 
than the onset of pa iring, T* (H) , in contrast to previous 
wor k) 3 | 6 | 7 | io | i4 | 25 | 26 | To tms en ^ ur work has established 

that the Gor'kov equations lend themselves to the non- 
linear, analytic approach required, provided only degen- 
erate energy states are paired. It has also shown that 
unique pairing partners are not required for a tractable 
theory. A robust result of this theory is that gapless 
states are present for both real and momentum-space 
pairing theories in a very high field. 

In summary, it is hoped that this formalism lays the 
foundation to explore a variety of magnetic field effects 
in the pseudogap phase and throughout the BCS-BEC 
crossover. 

We thank Victor Gurarie, Tin-Lun Ho, and Vivek 
Mishra for helpful discussions. This work is supported 
by NSF-MRSEC Grant 0820054. P.S. acknowledges sup- 
port from the Hertz Foundation. 



Appendix A: Details of the Orbit-Center Pairing 
Calculations 



Since v[/P air is independent of k z , that subscript will be 
dropped in pair terms that follow. 
Ref . [5] provides an identity, 

JVi+JV 2 

TpN u X+Y{ri)lpN 2 ,X-Y(r2) = ^ C p 



/ j ^N X N 2 
P=0 



X V ; p3s:( r cm)^5vJ+Ar 2 _P ! 2y( T *rel) 



(Al) 



where Ni,N 2 , and P are Landau levels, r cm = (r x + 
J*2)/2, r ro i = r\ — r-2, ^ cm is an orbital function with 
I — > 2~ 1 / 2 Z (appropriate for a charge-2e particle), ip rcl 
has I — ¥ 2 1 / 2 /, and (7jy N2 is a complicated combinatorial 
prefactor which is equal to I if Nx = N 2 = 0. Thus, we 
know that 



pair 



Ni +-/V2 
\ ' nP 



\,.\ 2 \,y( V ) ~ ^ ( ^N 1 N 2 ^p m x( r )' l l ; N i 1 +N2-P,2Y 



(0) 



p=o 



Nt+N 2 



L y L z 2 N ±+ N i- p {N 1 +N 2 - P)\ 



p=o v v 

=< (4)" v, """* S(r) 

We can now proceed to calculate the A mn elements for 
all possible pairs, with the result that 
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A> m =(N 1 ,x+Y), n =(N 2 ,x-Y) = \ dr A(r)^~^ fe , (r)^™°°y _ fcl (r) 



L x LyL z y/Trl 



- J dr^2 cxp (i^m^j V^ lbx (r) 



iV!+iV 2 / " . 1 . 1/4 



p=o 



A / s/2b x ~ t7r(b„/a)(X/b x ) 2 r / 1 f 1, p^/ig if X — rnh m (= J, 

A V L a L y L x ^iil a e L/ N 1 N 2 \J LyL ^ N i+ N 2{N 1+ N 2 )\ J e ltA - TOt^, TO € A, 

otherwise. 



Here <7° - , (-1)" 2 TO ^s™ 

nere Ojy lW - 2 - y jv^i 2(J v 1+ iv 2 )/2 to give 



[A / he , 1 e «7T (b y /a)(X b^y -V/ljj [fx-mb mf Z 

^TO=(iVi,X+K),n=(iV 2 ,X-r) — % v 

otherwise. 
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